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Outflows from the high-mass protostars NGC 7538 IRS1/2 
observed with bispectrum speckle interferometry^ 



Signatures of flow precession 
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r— s , Context. NGC 7538 IRS 1 is a high-mass (30 Mq) protostar with a CO outflow, an associated ultracompact H II region, and a linear methanol 
maser structure, which might trace a Keplerian-rotating circumstellar disk. The directions of the various associated axes are misaligned with 
each other. 

Aims. We investigate the near-infrared morphology of the source to clarify the relations among the various axes. 
Q Methods. A"-band bispectrum speckle interferometry was performed at two 6-meter-class telescopes — the BTA 6m telescope and the 6.5m 
V^ ' MMT Complementary IRAC images from the Spitzer Space Telescope Archive were used to relate the structures detected with the outflow at 
C/3 , larger scales. 

Cu ■ Results. High-dynamic range images show fan-shaped outflow structure in which we detect 18 stars and several blobs of diffuse emission. We 

^ I interpret the misalignment of various outflow axes in the context of a disk precession model, including numerical hydrodynamic simulations 

• "^ . of the molecular emission. The precession period is ~ 280 years and its half-opening angle is ~ 40°. A possible triggering mechanism is non- 

^^ ■ coplanar tidal interaction of an (undiscovered) close companion with the circumbinary protostellar disk. Our observations resolve the nearby 

, massive protostar NGC 7538 IRS2 as a close binary with separation of 195 mas. We find indications for shock interaction between the outflow 

■ activities in IRSl and IRS2. Finally, we find prominent sites of star formation at the interface between two bubble-like structures in NGC 7538, 

suggestive of a triggered star formation scenario. 

Conclusions. Indications of outflow precession have been discovered to date in a number of massive protostars, all with large precession 
angles (~ 20-45°). This might explain the difference between the outflow widths in low- and high-mass stars and add support to a common 
coUimation mechanism. 

Key words, stars: formation - stars: individual: NGC 7538 IRSl, NGC 7538 IRS2 - techniques: bispectrum speckle interferometry, interfer- 
ometric 

1 . Introduction cumstellar disks is still rare, whereas outflows seem to be 

omnipresent in the high-mass star forming regions. Outflows 
Protostellar disks and outflows are essential constituents of ^^^^^^ ^^^ ^^jy ^„g^j^^ momentum from the infalling mat- 
the star formation process. For high-mass protostellar objects ^^^^ ^^^ ^j^^ j^^jp ^^ overcome the radiation pressure limit to 
(HMPOs), direct evidence for the presence of compact cir- protostellar accretion, by carving out optica lly thin cavities 

along which the radiation pressure can escape (iKrumholz et alJ 

Send offprint requests to: skraus@mpifr-bonn.mpg.de l2005l) 

* Observations reported here were obtained at the BTA 6m tele- ., _ ,,■ , ■ ,- • , 

„,„.,. , ■ , r^. . , Ti,T>TT, How outflows are collimated is a matter or ongoing de- 

scope ot the Special Astrophysical Observatory and at the MMT ° ^ 

Observatory, a joint facility of the Smithsonian Institution and the ^ate and may depend on the stellar mass of the outflow- 
University of Arizona. In addition, this work is based in part on living source. One of tile arguments m support of this con- 
archival data obtained with the Spitzer Space Telescope, which is elusion is that outflows from high-mass stars appear less col- 
operated by the Jet Propulsion Laboratory, California Institute of limated than the outflow s and jets from their low-mass coun- 
Technology under a contract with NASA. terparts ( IWu et al.l2004l) . Therefore, it has been suggested that 
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outflows from HMPOs might be driven by strong steflar winds, 
lacking a recollimation mechanism. Since HMPOs typicaUy 
form in dense clusters, another possibility is confusion by the 
presence of multiple collimated outflows. 

However, since there is evidence that the binary frequency 
is significantly higher for h igh-mass than for low-mass stars 
(e.g., IPreibisch et al. 19991) . another possibility is that out- 
flows from HMPOs simply appear wider, assuming they un- 
dertake precession. A few cases where outflow precession 
have been proposed for HMPO o utflows (e.g. Shepherd et al. 
l200(i IWeigeltetal.ll2002l l2005l) show precession angles of 
~ 20 to 45°; considerably wider than the jet precession an- 
gles of ty pically just a few de grees observed towards low- 
mass stars ( iTerauem et al.l 19991) . This is in agreement with the 
general picture that high-mass stars form at high stellar den- 
sity sites and therefore experience strong tidal interaction from 
close companions and stellar encounters. 

The detection of precessing jet-driven outflows from 
HMPOs adds support to the hypothesis of a common forma- 
tion mechanism for outflows from low to high-mass stars. 
Furthermore, jet precession carries information about the 
accretion properties of the driving source and, simultaneously, 
about the kinematics and stellar population within its closest 
vicinity, yielding a unique insight into the crowded places 
where high-mass star formation occurs. 

In this paper, we report another potential case of outflow 
precession concerning the outflow from the high-mass (30 Mq, 
IPestalozzi et alJ2004 protostellar object NGC 7538 IRSl. 

We obtained bispectrum speckle interferometry of IRSl 
and IRS2, which provides us with the spatial resolution to 
study the inner parts of the outflow, detecting filigreed fine 
structure within the flow. Information about even smaller 
scales is provided by the intriguing methanol maser feature, 
which was detected at the position of this infrared source 
and which was model ed successfully as a p rotostellar disk 
in Keplerian rotation JPestalozzi et alJ 120041) . To search for 
outflow tracers on larger scales, we also present archival 
Spitzer/IRAC images. In addition, this allows us to relate the 
sources studied with bispectrum speckle interferometry with 
the overall star forming region and we find new hints for 
triggered star formation in this region. 



1.1. Previous studies of NGC 7538 

The NGC 7538 molecular cloud is located in the Cas 
OB2 associat ion in the Perseu s spiral arm at a distance of 
~ 2.8 kpc tel itz et alJ Il982l) . Several authors noted that 
NGC 7538 might present a case of triggered or induced star 
formation since it shows ongoing star formation at various 
evolutionary stages, apparently arranged in a northwest (most 
dev eloped) to southeast (youn gest evolutionary stage) gradi- 
ent (iMcCaughreaneTaillilll). 

At optical wavelengths, the appearance of the region is 
dominated by diffuse H II emission, which e xtends several ar- 
cminutes from the southeast to the northwest ( iLvnds & O'Neill 



Il986l) . In 1974, IWvnn-WiUiams. Beckhn. & Neugebaueil de- 
tected eleven infrared sources (IRS 1-11) in the NGC 7538 re- 
gion, wherein IRS 1-3 are located on the southeast-corner of 
the fan-shaped H II emission in a small cluster of OB-stars. 
IRS 1 is the brightest NIR source within this cluster and is em- 
bedded within an ultracompact (UC) HII region whose size 
was estimated to be ~ 0'.'4 ( rtp ^ 10^ cm~^ , measured in 
5 and 15 GHz CO continuu m ICampbellll 1 984 The spectral 
type was estimated to be 07 JAkabane & Kundl2005l) . which 
implies a luminosity ~ 9.6 x IQ'^Lq. VLA observations with 
a resolution down to 0*.'! (=180 AU) also revealed a double- 
peaked structure of ionized gas within the UC core (peaks sep- 
arated by ~ 0'.'2), which was in terpreted as a disk coUimating 
a nor th-south-oriented outflow ( ICampbellll 1984i iGaume et alJ 
I1995I) . This interpretation is also supported by the detec- 
tion of elongation of the d ust-emitting region at mid-infrared 
(MIR) wavelengths_(5jan: Hack wellet al.lll982t 1 1 .7 jum and 
18.3 urn: IPe Buizer & Minierii2005l) and imaging studies per- 
formed in the sub-millimete r continuum (350jum, 4 50 fim, 
800 ;um, 850 ;um, 1.3 mm: ISandeU & Sieversll2004 show- 
ing an elliptical source with a si ze of ~ 11"6 x 7" 6 along 
pa' ~ -80°) and CO line emission dScoville et al.l 19861 show- 
ing a disk-like structure extending ~ 22" in the east-west direc- 
tion). Also, polarization measurements of the infrared emission 
arou nd IRS 1 can be construed in favor of the d isk interpreta- 
tion (lpvck& Cappsll 19781 iTamuraet al]ll99ll) . iKawabe et ail 
( I1992I) carried out interferometric CS (J=2 ^1) observations 
and found a ring-like structure, which they interpret as a nearly 
face-on protostellar disk of dense molecular gas. 

Furth er evidence for outflow activity was found 
by Gau me et alJ ( 1 19951) . who measured the profile of the 
H66ff recombination line and derived high velocities of 
250 kms ', indicating a strong stellar outflow from IRSl. 
CO (J=l — > 0) spect ral line mapping showed a bipolar 
flow JFischer et alJl985l) . The mass outflow rate Mnnfflnw from 
IRSl was estimated to be ~ 5.4 x 10"^ Mp vr~' JDayis et all 
ll£98). Interferometric observations by IScoville et alJ 09861 
beam size 7") show that the blue and red-shifted lobes are 
separated by 28" with a position angle of -45°, and IRSl is 
located on this axis just between the lobes of this high-velocity 
(-76 to -37 km s ') CO outfl ow. In comparing the data ob- 
tained with various beam sizes ( ICampbelll984lKameva et alJ 
ll989), these seem to indicate a change in the position angle 
of the flow direction at different spatial scales, ranging from 
PA ~ 0° at a.'3, PA ~ -25° at 2", PA ~ -35° at 7", to PA -40° 
at 16". 

Within the immediate (~ 0'.'5) vicinity of IRSl, a 
large variety of masers has been discovered, including 
OH fPic keletal. 1982), H 2CO (formaldehyd e, iRots et all 
1981; Ho ffman etalJl2003h . NH3 (ammonia, iMadden et all 
ll986i) . CH3OH (methanol, five feafiir es A, B, C, D, E 
were detected at 6.7 an d 12.2 GHz: iMenten et alJ 19861: 
iMinier et all 1 19981 I2OO0I) '^NHj JJohnstonet alJll989l) . and 
H2O JKameva et alJ Il990l) . Some of the masers show only 
vague signs for a systematic alignment within linear ('^NH3, 
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from north to east. 
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PA 60°) or ring-like structures (H2O, methanol-maser fea- 
ture E). However, the methanol-maser feature A represents one 
of the most convincing cases of systematic alignment, in both 
linear spatial arrangement (PA ~ -62°) and well-defined veloc- 
ity gradient, observed to date in any maser source. The qualita- 
tive interpretation of this structure as an edge-on circumstellar 
disk ( Minier et al 119981) was later c onfirmed by the detailed 
modeling of IPestalozzi et alJ 1120041) . which showed that the 
alignment in the position-line-of-sight (LOS) velocity diagram 
of maser feature A can be modeled accurately assuming a pro- 
tostellar disk with Keplerian rotation. 

Aiming for a more complete picture, seve ral authors (e.g. 
iMinier et alii 1 9981: lOe Buizer & Minieill2005l) also tried to in- 
corporate the presence of methanol maser features B, C, D, 
and E in the circumstellar disk model for feature A and inter- 
preted them as part of an outflow which is oriented perpedicu- 
lar to feature A. Since these maser features are southwards of 
the putative circumstellar disk, it remains unclear whv thev ap - 
pear blue-shifted with respect to feature A JMinier et alJl998l) . 
whereas the southern lobe of the CO-outflow is red-shifted. 

Besides the circumstellar disk interpretation for the origin 
of the maser featu re A mentioned above, an alt ernative scenario 
was proposed by IPe Buizer & Minien ( 120051) . who suggested 
that feature A might trace the walls of an outflow cavity. 

The region was also intensively observed in the infrared. 
Survey i mages of the infrared continu um emissio n were pre- 
sented bv lCampbell & PerssonI (1 19881 H, K) and lOjha et all 
1I2OO4 J, H, Ks) and showed diffuse emission, which ex- 
tends from the IRS 1-3 cluster in a fan-shaped structure to- 
wards the northeast and north, approximately tracing the op- 
tical H II region. The northeast border of this NIR emitting re- 
gion also appears very pr onounced in t he con tinuum- subtracted 
H2 2.122 jum maps bv IPavis et al.l lll998l) . possibly tracing 
the illuminated surfaces of nearby molecular clou ds or the in- 
ner w alls of a vast outflow cavity. Furthermore, IPavis et alJ 
( 119981) discovered two bowshock-shaped structures, centered 
roughly on the IRS 1-3 cluster and orientated again along the 
northwest-southeast direction (PA ~ -30°) in H2 2.122 //m. 
With imaging at arcsecond resolution and the use of sev- 
eral spectral filters (J, H, K, [Fe 11] 1 .65 urn. Bry 2.165 //m, 
H2 2.122 yum, and 3.29 um). iBloomer et al. (1998) attempted 
to identify the source and mechanism of the outflow. Based 
on a cometary-shaped morphology in the [Fe II] line images 
and shell-like rings observed in the J, H, and A'-bands, these 
authors propose a stellar wind bowshock model in which the 
motion of IRS2 relative to the molecular cloud produces the 
diffuse NIR emission within the vicinity of the IRS 1-3 cluster 

The first /T-band speckle images, ta ken with the 3.5 m- 
telescope on Calar Alto were presented bv lAlvarez et alJ (120041) 
and showed substructure in the vicinity of IRS 1 ; namely, two 
strong blobs (A, PA ~ -45°; B, PA ~ -70°), a diffuse emission 
feature (C, PA ~ 0°) as well as several faint point-like sources 
(«-/). 



2. Observations 

2. 1. Bispectrum speckle interferometry 

The first set of observations was performed on 2002-09-24 
using the 6.0 m BTA (Big Telescope Alt-azimuthal) tele- 
scope of the Special Astrophysical Observatory located on 
Mt. Pastukhov in Russia. Additional data were gathered 2004- 
12-20 with the MMT (Multiple Mirror Telescope) on Mt. 
Hopkins in Arizona, which harbors a 6.5 m primary mirror. 
As detector, we used at both telescopes one 512x512 pixel 
quadrant of the Rockwell HAWAII array in our speckle cam- 
era. All observations were carried out using a TT'-band fil- 
ter centered on the wavelength 2.12 //m with a bandwidth 
of 0.21 fim. During the BTA observation run, we recorded 
420 speckle interferograms on NGC 7538 IRSl and 400 in- 
terferograms on the unresolved star BSD 19-901 in order to 
compensate for the atmospheric speckle transfer function. The 
speckle interferograms of both objects were taken with an ex- 
posure time of 360 ms per frame. For the MMT observa- 
tions, the star 2MASS 23 134580H-6 124049 was used for the 
calibration and 120 (200) frames were recorded on the target 
(calibrator) with an 800 ms exposure time. The modulus of 
the Fourier transform of the object (vis ibility) was obt ained 
with the speckle interferometry method dLabevridl 197(1) . For 
image reconstructi on we used the bispectrum speckle inter- 
ferometry method (IWeigeltlll977l Iweigelt & Wirnitzeilll983l 
iLohmann et ai]ll983[ iHofmann & Weigeltlll986l) . With pixel 
sizes of 27.0 mas (BTA) and 28.7 mas (MMT) on the sky, the 
reconstructed images possess fields of views of 13'.'8 (BTA) 
and 13'.'1 (MMT), respectively. 

We found that the BTA data allows the highest spatial res- 
olution (and is therefore perfectly suited for the identification 
of point-sources within the field), whereas the image recon- 
structed from the MMT data allows a high dynamic range in 
the diffuse emission. Therefore, we show the diffuse emission 
within an image of moderate resolution (reconstructed from 
MMT data, see Figure Q];) and perform point-source identifi- 
cations within the higher resolution image reconstructed from 
BTA data (Figure^). In order to distinguish point-sources and 
diffuse structures reliably, we reconstructed images of various 
resolutions (146 mas, 97 mas, 72 mas) and carefully exam- 
ined changes in the peak brightness of the detected features. 
Whereas for point-sources the peak brightness increases sys- 
tematically, it stays constant or decreases for diffuse structures. 

To perform an absolute calibration of the astrometry in our 
images, we measured the position of IRSl and IRS2 in the Two 
Micron All Sky Survey (2MASS) Ks Adas images and use the 
determined absolute positions as reference for our astrometry. 
We estimate that the accuracy reached in the relative astrometry 
is ~ O'.'l. The absolute calibration introduces further errors (~ 

a:2). 



2.2. Spitzer/IRAC Archive data 

In order to relate our high-resolution images with the mor- 
phology of the NGC 7538 molecular cloud at large scales, 
we examined archival 3.6, 4.5, 5.8, and 8.0 jum images (PI: 
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Fig. 1. Bispectrum speckle images 
(/T'-band) reconstructed from data 
taken with a) the 6.5m MMT and 
b) the 6m BTA telescope. To show 
the weak emission features, the in- 
tensity of IRSl was clipped to 
2% of the total flux. Within the 
high-resolution image {b), speckle- 
noise artifacts appear around IRSl 
(marked with a circle). These weak 
features represent small distortions 
of the point-spread-function (PSF) 
on the 1%-level and do not influ- 
ence the reliability of the identifi- 
cation of point sources within the 
image. The absolute coordinates of 
IRSl are a = 23*" 13" 45^35 and 
6 = -h61° 28' 10'.'84 (J2000, de- 
termined from 2MASS, accuracy ~ 
a.'S). 
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G. G. Fazio), tak en with the Infrared Array Camera (IRAC, 
iFazio et al J 120041) on the Spitzer Space Telescope. The four 
bands are recorded simultaneously using two InSb (3.6 jum, 
4.5 jum) and two Si: As (5.8 yum, 8.0 yum) detectors. The central 
wavel engths and bandwidths of the IRAC bands ( iHora et alJ 
l2004 are 3.56 //m (A^ = 0.75 //m), 4.52 //m (A^ = 1.01 //m), 
5.73 jum {^A = 1.42 ^im), and 7.91 ^lm {^A = 2.93 //m). 
Each image consists of 256 x 256 pixels, corresponding to 
a ~ 5' X 5' field-of-view on the sky. The data used include 
48 Spitzer pointings taken on 2003 December 23 in the High 
Dynamic Range (HDR) mode. In HDR mode, for each point- 
ing, images are taken with two exposure times (0.6 s and 12 s) 
in order to record both bright and faint structures. However, 
the two brightest sources, IRSl and IRS9, are saturated even 
within the 0.6 s exposure. 

We used the mopex software (2005-09-05 version), released 
by the Spitzer Science Center (SSC), to process both the long 
and short exposure images. Beside the basic calibration steps 
appHed by the Basic Calibrated Data (BCD) pipeline (S 1 1 .0.2), 
we performed Radhit detection, artifact masking, and pointing 
refinement. Finally we generated a mosaic in which the satu- 
rated pixels of the long exposure image were replaced by the 
corresponding pixels of the 0.6 s exposure. The optical design 
of IRAC induces a shift of ~ 6.'8 between the 3.6/5.8 fim and 
4.5/8.0 jim pointings, leaving an overlap of 5 ! 1 between all four 
bands. 

In Figure 121 color composites of the 3.6/4.5/5.8 /im and 
4.5/5.8/8.0 yum band images are shown. 

The diffuse emission in three of the four IRAC bands 
is dominated by Polycyclic Aromatic Hydrocarbons (PAHs, 
IChurchwell et al. 2004), which trace the border of regions 
excited by the UV photons from HMPOs particularly well. 
Contributions are also expected from several vibrational lev- 
els of H2 (Smith &_Rosen 2005b), atomic lines, CO vibrational 
bands, and thermal dust grain emission. 



3. Results 

3. 1. Bispectrum speckle interferometry: Small-scale 
structures around IRS1/2 

3.1 .1 . IRS1 Airy disk elongation and diffuse 
emission 

In our speckle images, the Airy disk of IRSl itself appears 
asymmetric, being more extended towards the northwest di- 
rection (PA ~ -70°, see Figure^ and inset in the lower left 

of Figure|3}- In the same direction (PA 60°), we find two 

strong blobs (A, B + B') of diffuse emission at separations of 
~ 1" and 2". These blobs and additional diffuse emission seem 
to form a conical (fan-shaped) region with a 90° opening angle 
extending from IRS 1 towards the northwest. 

A careful examination of the power spectrum of IRSl has 
shown that the detected asymmetry of IRSl is not caused by a 
companion, but seems to represent diffuse emission. Therefore, 
we can rule out a close binary system of similar-brightness 
components down to the diffraction limit of ~ 70 mas. For the 
case of a binary system with components of significantly dif- 
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Fig. 4. By measuring the speckle noise around the PSF of IRS 1 , 
we can rule out binarity of IRS 1 on a 3cr-level as a function of 
apparent separation and intensity ratio. 

ferent brightnesses, we can put upper limits on the brightness 
of the hypothetical companion as a function of the projected 
binary separation (see Figurel^. 

The PA of the elongation of the Airy disk is similar to the 
PA of K' blobs A, B, and B' . Another strong feature (C) can 
be seen towards PA~ 0°. The blobs seem to be connected by 
a bridge of diffuse emission extending from feature B to C. 
Overall, the diffuse emission seems to form a fan-shaped re- 
gion which is extending from IRS 1 towards the northwest with 
an opening angle of nearly 90° . We identified some further fea- 
tures and list their position angles and separations in Table ^ 
The directions, which were reported for various outflow trac- 
ers, are also listed in this table and illustrated in Figure|5] 

Our features A, B + B', and C appear to coincide with the 
features A, B, and C identified bv lAlvarez et al..t.2004.) . A com- 
parison suggests that features A + B + B\ C. D cor respond to 
the peaks IW, lA^, and 1 A^ W in Ffamura et al J (1 1 99 ih . 

3.1 .2. Binarity of NGC 7538 IRS2 

IRS2 appears resolved as a close binary system. Using an im- 
age reconstructed from BTA data with a spatial resolution of 
80 mas (see inset in the upper-left corner of Figure |3}, we de- 
termined the separation to be 195 mas and found a PA of -123° 
for the 1™9 fainter companion (2002-09-24). We designate the 
brighter component in the /T'-band as IRS2fl and the fainter as 
lRS2b. 



3.1 .3. Detection of fainter cluster members 

Besides IRSl and IRS2fl/i', we were able to identify 18 ad- 
ditional fainter point-like sources (a-r) within the BTA image, 
whose positions and /T'-band magnitudes are listed in Table|2l 
In order to test whether these sources are physically re- 
lated to NGC 7538, one can compare the stellar number density 
for the brightness range 11™0 to 12'."0 in our speckle image, 
A''speckie = 18/128 arcsec^ a; 2.1 x 10''/deg^, with the num- 
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Table 1. NGC 7538 IRSl outflow directions reported for various tracers. 
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Structure 


Beam Size 


Scale 


PA of 


Dynamical 


Ref. 


Comments 








["] 


["] 


outflow 
direction [°] 


Age^' 
[10^^ yrs] 




(Velocities in [km s"']) 




Methanol masers 




0.03" 




+ I9M 


- 


[8] 


see Fig. 13 




1.0 cm cont. 


inner core < 0'.'5 


0.13 


0.4 


-^0" 


<0.03 


[9] 


/, = 0.15 yrs" 




1.3 cm cont. 


inner core < 0'.'4 


0.11 


0.4 


-^0" 


<0.02 


[5] 






6.0 cm cont. 


inner core < 0'.'75 


0.35 


0.4 


-1-0 


<0.04 


[1] 






1 .0 cm cont. 


outer core > 0'.'5 


0.13 


1.0 


-25" 


>0.03 


[9] 


t^ = 0.3 yrs" 




1.3 cm cont. 


outer core > 0'.'4 


0.11 


1.0 


-25" 


>0.02 


[5] 






6.0 cm cont. 


outer core > 0'.'75 


0.35 


1.0 


-15... -20 


>0.04 


[1] 






6.0 cm cont. 


outer core > 0'.'5 


0.4 


1.0 


-25" 


>0.03 


[9] 






MIR11.7;[/m 


IRS 1 Elongation 


0.43 


3 


-45 


0.16 


[10] 


see Fig.|6}: 




MIR 18.3y[/m 


IRS 1 Elongation 


0.54 


4 


-45 


0.16 


[10] 


see Fig.|6j 




NIRAT'-band 


IRS 1 Elongation 


0.3 


0.6 


-78 


0.03 


- 


see Fig.|3l 




NIRif'-band 


feature A 


0.3 


1.6 


-65 


0.08 


- 


see Fig.|3l 




NIR K'-hmd 


feature F 


0.3 


2.1 


-33 


0.11 


- 


see Fig.|3l 




NIR a: '-band 


feature B 


0.3 


3.0 


-39 


0.16 


- 


see Fig.|3l 




NIR a: '-band 


feature B' 


0.3 


3.3 


-57 


0.18 


- 


see Fig.|3l 




NIR a: '-band 


feature C 


0.3 


4.8 


-1-6 


0.25 


- 


see Fig.|3l 




NIR a: '-band 


feature E 


0.3 


6.2 


-20 


0.33 


- 


see Fig.|3l 




NIR a: '-band 


feature D 


0.3 


7.4 


-HlO 


0.39 


- 


see Fig.|3l 




NIR a: '-band 


eastern wall 


0.3 


- 


-H25 


- 


- 


see Fig.|3l 




NIR a: '-band 


western wall 


0.3 


- 


-65 


- 


- 


see Fig.|3l 




[Fell] 1.65 yum 




1 


15° 


N-S" 


0.8 


[7] 


around IRS2; see Fig.|6j 




H2 


northern bow 




30° 


N-S 
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see Fig.|6|f 
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Note - V,- and Vi, are measured relative to the velocity of methanol maser feature A ( V = V - 56.25 km s" ) 

" Estimated from figures presented within the reference paper; therefore, with limited accuracy. 

'' The half-separation between the red- and blue-shifted CO lobe is given. 

" For VLBI observations, we give the estimated error on the absolute position of the measured maser spots. 

'' The expected outflow direction is given; i.e., perpendicular to the measured orientation of maser feature A. 

^ Electron recombination time given in the reference paper. ^^^_^^^^ ^^^ 

^ Assuming an outflow velocity of 250 km s ' , which was measured bv lCaume et alj J 19951) within the 11660- recombination line. For the 
CO emission, we also use the measured CO outflow velocity and provide the corresponding dynamical age in brackets. Since all velocities are 
measured alon e LOS, this timesc ale can only provide upp er li mits. ^_^ ^^^^^^^^ ^_^ ^^^_^^^^ 

Re ferences: [l] |Campbell"(1984); [2] Fischer et al. ( 1985); [3] ScoviOe et al. ( 1986); [4] Kameva et al. ( 1989); [5]"G aume et"ai]jl995h ; 
[6] lDavis et a'illl99S!); [71 Bloomer et alj lil998) ; [81 Minier et al. (2000) ; [91 Franco-Hernandez & Rodriguez l2004i) ; 



[lOl lDe Buizer & Minieij i2005h 



71J 



ber expected from the cumulative /T-band luminosity function 
(KLF)ofthe NGC 7538 field- A^fieid ~ 1.8xlOVdegl Although 
these number densities were obtained with different spatial res- 
olution, the clear over-density of stars in our speckle image is 
significant and we conclude that most of the detected stars are 
likely members of the NGC 7538 star forming region. When 
using the KLF for the IRS 1-3 region instead of the whole 
NGC 7538 field, the stellar over-density in our speckle image 



^ The A'-band luminosity function by iBalog et aO yOOJ) for the 
whole NGC 7538 region, corrected with the on-cluster KLF, and cu- 
mulated for the magnitude range 1 1'^O to 12'^0 was used. 



becomes even more evident (Mrs 1-3 cluster ~ 1-4 x 10^/deg^). 
Since these stars are about 5 to 6 magnitudes fainter than IRSl, 
they are likely to be part of the associated intermediate mass 
stellar population. 

The arrangement of the stars within the fan-shaped nebula 
does not appear to be random, but follows the S-structure of the 
diffuse blobs (see Figure 1^. Most remarkable, more than half 
of the stars seem to be aligned in a chain reaching from feature 
B to C (PA ~ 45°). Within the diffuse blobs close to IRSl (A, 
B, B'), no stars were found, whereas embedded in blob C, three 
stars could be detected. 
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Table 2. Point sources identified in our speckle images. For details about IRSl and the binary system lRS 2a/b, we refer to the 
text. We identify components a to / with the stars already discovered in the image by JAlvarez et alJ ( 120041) . 



Name RA (J2000y' DEC (12000)" A:' Magnitude'' Comment 



a 


23" 13" 45'!33 


61° 28' 21'.'02 


10^94 




b 


44^:81 


ia'63 


1 1?29 




c 


44^:95 


12'.'33 


11™46 




d 


45U1 


09'.' 69 


11"^73 




e 


45^37 


15'.' 12 


11-^77 


embedded in feature C 


f 


4S34 


15'.'93 


ir73 


embedded in feature C 


g 


4S36 


14'.' 80 


ll-?73 


embedded in feature C 


h 


4S31 


14'.'72 


11?86 




i 


45'=.29 


13'.' 74 


ir71 




J 


4S27 


13'.' 25 


ir73 




k 


4S25 


13'.' 05 


ir8i 




1 


45^:24 


13'.' 28 


ir83 




m 


45^:22 


12'.'71 


11?71 




n 


4S36 


12'.' 99 


11?75 







45^:50 


10'.'50 


11"^ 60 




P 


4S29 


18'.' 42 


11?73 




q 


4y;i6 


19".' 69 


11?77 




r 


45^33 


08'.' 48 


11-?71 





" For the astrometry, the relative errors are of the order of 0'.' 1 . The absolute calibration using the reference position of IRS 1 in 2MASS 
introduces further errors (0'.'2). 

* The photom etry was done rela tive to IRS 1 with an uncertainty of 0"3. For the conversion to absolute photometry, we assumed a IRS 1 
magnitude of 8™9 JOiha et all2004l) . 



3.2. Spitzer/IRAC: Morphology at large spatial 
scales 

Imaging of NGC 7538 at optical wavelengths showed that 
diffus e emission c an mainly be found in the vicinity of 
IRS5 jLynds & O'Neil 198^1) . At near-infrared (NIR) wave- 
lengths JOiha et alJl2004l) . a diffuse structure can be found ex- 
tending from the IRS 1-3 cluster towards the northwest with the 
strongest emission around IRS5. 

The Spitzer/IRAC images reveal a more complex, bubble- 
like structure (see Figure ^, c), whose western border is 
formed by a pronounced ridge-like filament connecting IRSl- 
3 with IRS4 and reaching up to IRS5 (see Figure |2j;). At the 
western border of the bubble a wide conical structure is lo- 
cated, with a vertex on 2MASS 23 135808H-6 130484. Another 
conical structure can be detected close to the northern bor- 
der of the bubble. Several other outflow structures can be 
found in the IRAC image; most noteworthy, the unidirec- 
tional reflection nebula around 2MASS23144651H-6129397, 
2MASS2313169lH-6129076,and2MASS23130929H-6128184 
(see FigureEt). The sources 2MASS23131660H-6128017 and 
2MASS 23133184H-6125161 appear to be embedded in a shefl- 
like cloud. 

Besides the position of the strongest near-infrared sources. 
Figure |2] shows also the position of the submi l limete r 
(450 yum, 850 yum) clumps reported bv lReid & Wilsonl ( l2005l) . 
These clumps trace the filaments and knots of the bubble, 
which can be seen in the IRAC images, very well. Besides 
this, the submillimeter clumps suggest another bubble-like 
str ucture to the southw est of IRS4 (see also the images 
in iReid & Wilsonl 120051) . This bubble seems to be invisible 



at near- and mid-infrared wavelengths, although several NIR 
sources are located on its border (2MASS 23130929H-6128184, 

2MASS23133184h-6125161).^ 

As already pointed out by iReid & Wilsonl ( l2005l) . it is in- 
teresting to compare the position of the detected H2O masers 
with the position of the centers of high-mass star formation in 
the region and to find agreement in many cases (IRS 1-3, IRS9, 
NGC 7538S). However, as can be seen in Figure |2J for four 
locations of H2O masers, no MIR counterpart can be found in 
the IRAC images (the detection limits for point sources in the 
four IRAC bands are roughly 3.6, 5.3, 31, and 34 fi}y for the 
IRAC bands at 3.6, 4.5, 5.8, and 8 jum assuming medium sky 
background). 

4. Discussion 

4. 1. Nature of the observed K'-band emission 

In the wavelength range of the ^'-band filter (Aq = 2.12 fim, 
AA = 0.21 //m), we record not only continuum radiation 
(e.g. scattered light, thermal dust emission, stellar contin- 
uum emission), but also l ine emission (e.g. H? ). However, 
both [ Bloomer et alJ lll998L see Figure |6jl and IPavis et alJ 
II1998I) did not detect significant amounts of H2 emis- 
sion around IRSl. F urt hermore, deep optical ima ging by 
lElsaesser et"ai] ( Il982l) and lCampbell & Perssoi] (Il988l) reveal a 
weak optical source offset 2'.' 2 north of the radio source IRSl. 
The latter authors argue that the strong extinction (Ay = 13) 
derived for IRS 1 makes it highly unlikely that this optical emis- 
sion is connected to IRS 1 itself but that it most likely represents 
scattered light. The measured offset suggests that the faint op- 
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tical source should be associated with blobs A and B in our im- 
ages, making scattering the most likely radiation mechanism 
for the detected /iT'-band emission. This conclusion is also sup- 
ported by polarization measurements (Dvck & Lonsdale 1979) 
which show a strong polarization of the 2 //m emission, trac- 
ing either scattered light or light transmitted through aligned 
grains. Henceforth we presume continuum to be the most im- 
portant contributor to the detected emission. 

4.2. Methanol maser feature A: Protostellar disk or 
outflow? 

We note that the 2MASS position of IRSl 
(a = 23M3™ 45^.35, 5 = +61° 28' KT-'S, J2000) and the po- 
sition of the methanol maser feature A (a = 23^ 13 "' 45^.364, 
6 = +61° 28' l(y.'55. .12000) reported bv lMinier et al. (2000) 
coincide within the errors^. Therefore, the methanol masers 
and the outflow driving source are likely causually connected, 
however a random coincidental alignment cannot be ruled out. 

Since methanol masers can trace both protostellar disks 
and outflows, it is not a priori clear how the linear alignment 
of the methanol maser feature A and the observed veloc- 
ity gradient should be in terprete d. For IRS l, both claims 
have been made (Pestalozzi et al.l 12004: De Buizer & Minien 
l2005l) . However, detailed modeling has provided strong 
quantitative support for the disk interpretation but is still 
missing for the outflow interpretation. Furthermore, a study 
by Pestallozi et al. (in prep.) suggests that simple outflow 
geometries cannot explain the observed properties of feature A. 

A major difference between these two scenarios is the 
orientation of the disk associated with the outflow driving 
source: Whereas in the disk scenario the methanol masers 
are lined up within the disk plane (PA ~ -62°), the outflow 
scenario suggests an orientation of the disk plane perpendicular 
to the maser alignment (PA ~ +28°). The observed asym- 
metry in our NIR speckle images, as well as the elongation 
of the emission observed in the 11.7 and 18.3 yum images 
by IPe Buizer & Minier (2005i) . can be explained within both 
scenarios: 

Scenario A: If maser feature A traces an outflow cavity, the 

detected asymmetry might simply reflect the innermost walls 
of this cavity (oriented northwest), whereas the southeastern 
cavity of a presumably bipolar outflow might be hidden due to 
inclination effects. 

Scenario B: Alternatively, if the masers trace an edge-on 
circumsteUar disk, the asymmetry of the infrared emission 
could trace the western wall of an outflow cavity with a wide 
half-opening angle. The asymmetry cannot be attributed to the 
disk itself because the detection of stellar radiation scattered 
off the disk surface at such a large distance is highly unlikely. 



For completeness, we also mention the interpretation 
by iKameva et al.l ( Il989l) . who attributed the change between 
the direction observed in the UCHII region (PA 0°) and the 
high-velocity CO flow (PA -60°) to flow deflection, either by 
large-scale magnetic fields or due to density gradients. 

We proceed now to discuss both scenarios within an 
outflow-cavity model (Sec. 14. 3> and a precessing jet model 
(Sec. I4.4> . incorporating the large amounts of evidence col- 
lected by various authors over the last three decades. 

4.3. Scenario A: Outflow cavity model 

Since the intensity of the diffuse emission in our images seems 
to decrease with distance from IRS 1 and the vertex of the fan- 
shaped region appears centered on IRSl, we cannot support 
the hypothesis by Bloomer et al. ( 1998) , who identified IRS2 
as the likely source of the diffuse NIR emission. Instead, the 
observed fan-shaped region can be interpreted as a cavity that 
was formed by outflow activity from IRS 1 . Because the walls 
of the fan-shaped structure are well-defined, we can measure 
the opening angle of the proposed outflow cavity from the east- 
ern wall (PA 25°) to blob A (PA -65°), obtaining a wide total 
opening angle of 90°. 

The unidirectional asymmetry of IRS 1 in the NIR and MIR 
images (see Figure |6j is naturally explained in this context as 
scattered hght from the inner (< 1500 AU) walls of the cavity. 
This scenario is also consistent with the southeast-northwest 
orientation of the CO outflow, aligned roughly parallel to the 
methanol masers (PA ~ -62°). Blobs A, B, B' are located 
within the same direction and might resemble either clumps in 
the cavity or recent ejecta from the outflow. The various blobs 
might also indicate the presence of several outflows. 

In order to resolve the misalignment of the radio-continuum 
core with respect to the other outflow tracers, it was proposed 
that the radio-contin uum emission mig ht arise from a photoe- 
vaporated disk wind (iLugo et al.l2004l) . 

However, as noted above, the methanol maser feature A 
lacks a quantitative modeUng up to now. 

4.4. Scenario B: Precessing jet model 

A.A.^ . Constraints from the methanol maser disk 



^ The astrometric accuracy of the 2MASS 

catalogue was reported to be ~ Ol'lS (see 

http://ipac. caltech.edu/2mas s/releases/allsky/doc/expl-sup.htmT[i. 



The ci rcumstellar-disk modeling presented bv lPestalozzi et al.l 
(120041) reproduces the observational data for maser feature A in 
minute detail. Assuming a central mass of 30Mo and Keplerian 
rotation, this model confines the inner (r, ~ 290 AU) and 
outer {ro ~ 750 AU) radii of the disk (these radii scale as 
(M/30Mo)"'^^ with the central mass M). The model does not 
set strict constraints on the inclination and orientation of the 
disk on the sky. 

An uncertainty in the disk inclination arises from the as- 
sumption that methanol is formed within a surface layer of the 
disk from photoevaporation of H2O. The midplane of the disk 
might therefore be inclined within certain limits. A distinct in- 
clination is suggested by the fact that the NIR/MIR continuum 
emission, as well as the H2 shock tracer line emission, appears 
more pronounced towards the north than towards to south (see 
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Figures n and |5Ji). An inclination of the northern outflow to- 
wards us is also ind icated by CO outflow observations (e.g. 
iKameva et alJI 19891) . which show the blueshifted CO lobe of 
IRSl towards the northwest (see Figures|6lr &|6t). 

The disk orientation on the sky can only be constrained by 
the maser observations with a limited accuracy since the masers 
only trace a narrow latitudinal arc of the disk, missing potential 
disk warping. Nevertheless, it is still reasonable to identify the 
disk orientation with the linearly-aligned feature A. 

4.4.2. Indications for disk and jet precession 

Assuming that the alignment of the masers is representative of 
the orientation of the disk midplane (i.e. assuming disk warp- 
ing is negligible"*), it is evident that the direction perpendic- 
ular to the disk plane (PA -1-19°, the expected outflow direc- 
tion) is significantly misaligned from the axes of the bipolar 
CO outflow and the NIR fan-shaped region (PA ~ - 20°; illus- 
trated within Figure E). A lso, the observed bending dCampbelll 
Il984i lOaume et alJll995l) in the radio continuum could indi- 
cate a change in the outflow direction. Whereas the inner core 
(< (y.'5) is orientated along PA ~ 0°, the outer core (U.'S - 
r.'O) bends sUghdy towards the west (PA ~ -25°). This might 
indicate that the outflow changed its direction by this amount 
within the times needed by the jet to propagate the appropriate 
projected distances (~ 25 and ~ 50 years). 

The bending detectable in the UC H 11 region on scales of 
< I'.'O seems to continue at larger scales within the morphology 
observed in our speckle images, suggesting an S-shaped fine- 
structure of the diffuse emission extending from IRS 1 initially 
towards the northwest and further out towards north. The blobs 
A, B, and B' observed close to IRSl (PA 60°) might rep- 
resent the most recent ejecta, whereas the weak features which 
appear further away in our images (C, D) might trace earlier 
epochs of the history of the outflow. 

Based on these indications, we suggest a disk and jet pre- 
cession model. The fan-shaped diffuse emission in which the 
S-structure is embedded can be explained as scattered light 
from the walls of an outflow cavity, which was cleared by the 
proposed wandering jet. 

The western wall of this wide, carved-out outflow cavity 
might appear within our NIR and the MIR images as an elon- 
gation of IRS 1 . Since this elongation extends mainly towards 
the northwest, there must be an additional reason why the west- 
ern wall of this cavity appears more prominent than the east- 
ern wall. A possible explanation might be shock excitation of 
the western wall, which would cool through emission in shock 
tracer lines like H2, which is contributing to the recorded K'- 
band. 

Assuming the precession period derived in Sec. 14.4.31 the 
outflow (which currently points towards PA ~ 19°) would have 
excited the western wall of this cavity ~ 140 years ago, which 
corresponds roughly to the H2 radiative cooling time. 



The arrangement of the fainter cluster members embedded 
within the diffuse emission can be understood in this context, 
too: Taking into account that IRS 1 is still deeply embedded in 
its natal circumstellar cloud, the jet would have cleared the en- 
velope along its wandering path. The decreasing column den- 
sity results in lower extinction along the jet's path, revealing 
the fainter stars which likely formed in the vicinity of IRS 1 . 
The fainter stars might therefore be detectable only in those 
regions where the precessing jet reduces the extinction suffi- 
ciently. Within the blobs closest to IRSl, stars may be unde- 
tectable because of either inclination effects or confusion with 
the significantly higher surface brightness of blobs IRS 1 A, B, 
and B' (limiting the sensitivity to detect point sources), or be- 
cause of the high density of the outflowing material itself, pro- 
viding intrinsic extinction. 

The outflow tracers observed at rather large scales (CO, 
H2, see Figure |6j are oriented roughly in the same direc- 
tio n as the NIR f an-shaped structure. The CO channel maps 
by IScoville et alJ 111986. Figure |6^) suggest a change in the 
orientation of the CO outflow lobes for low and high veloci- 
ties. Whereas the low velocity CO outflow is oriented along the 
east-west direction, the high velocity lobes are oriented along 
PA ~ -35°. As CO traces material swept-up by the outflow and 
has a relatively long cooling time (of the order of lO"* yrs), the 
different orientations observed at low- and high velocities are 
more difficult to interpret. 

Finally, we speculate that the precession model might also 
explain why the velocities of the methanol maser features B, C, 
D, an d E are in the same range a s the velocities of the CO out- 
flow JDe Buizer & Minieiil2005l) . but show opposite signs for 
the LOS velocity with respect to feature A (maser features B, 
C, D, and E are blue-shifted, whereas the southern CO lobe is 
red-shifted). Assuming precession, the CO outflow would trace 
the average outflow direction around the precession axis (with 
the southern axis oriented away from the observer), whereas the 
methanol masers might trace clumps very close to the source, 
which were excited more recently when the southern part of the 
outflow was pointing towards the observer^. 

In general, precession can explain the change in the flow 
orientation, but potential alternative explanations include den- 
sity gradients in the surroundings of IRSl, the presence of mul- 
tiple outflows, and flow deflection. 

4.4.3. Analytic precession model 

In order to get a rough estimate for the precession parameters, 
we employ a simple analytical model with constant radial out- 
flow speed V. On the radial motion we superpose a precession 
with period Pprsc leading to the wave number 

2n 



k = 



vP, 



(1) 



prec 



"* Interestingly, the converse assumption (that disk warping is non- 
negligible) implies disk precession as well, as the jet would be 
launched perpendicular to the warped surface of the inner part of the 
disk. 



by the time that ejected material travels a distance r from the 
source, the direction of the jet axis changes by the angle kr. 

' This is consistent with the precession parameters determined in 
Sect. l4.4.3l where we find that the half-opening angle of the precession 
is larger than the inclination of the precession axis with respect to 
the plane of the sky if/. 
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To describe the jet propagation in three dimensions we in- 
troduce a Cartesian coordinate system centered on IRS 1 whose 
z-axis is along the line of sight (see Figure^ top). The preces- 
sion axis is in the y-z plane inclined by angle tfr to the plane of 
the sky^, and the jet axis makes an angle with it. For counter- 
clockwise^ precession, the coordinates of material at distance r 
from the origin are 



- rx 



^ sin© cos 
cos cos i/* + sin sin (p sin ifr 
- cos © sin i/* + sin sin (p cos tp ^ 



(2) 



where </> = (po + kr is the jet's azimuthal angle from the x- 
axis. The initial phase (pQ can be taken as 0° since the direction 
perpendicular to the methanol feature A seems to coincide with 
the eastern wall of the outflow cavity. The PA a of the y-axis 
can be set from the average angle of the fan-shaped region in 
the speckle i mages as a = -30 °. Using v - 250 km s"', as 
measured bv lOaume et al.l 1119951) within the H66a recombina- 
tion line, leaves as free parameters f prec, ©, ^, and the sense 
of rotation. Trying to fit the orientation of the maser disk, the 
orientation of the UC H 11 region, and the position of the NIR 
blobs with these parameters simultaneously, we find reason- 
able agreement with a precession period Ppi-ec = 280 ±10 yrs, a 
precession angle © = 40° ± 3°, a counter-clockwise sense of ro- 
tation, and small inclination i^ = 0° + 10°. At larger inclination 
angles (> 10°) loops start to appear, significantly degrading the 
agreement. In Figure^ we show the projected trajectory of the 
proposed wandering jet with the thick line, whereas the thin 
lines give the path obtained with a variation of +5° in a and if/, 
resembling the finite width of the flow. 

The analytic model presented in this section might suffice 
in order to get a rough estimate of the precession parameters, 
although it does not take into account the interaction of the flow 
with the ambient medium nor the excitation and cooling of the 
ambient material. 

These parameters can be used to predict how the orienta- 
tion of the methanol maser disk changes with time. Using the 
PA at the phase 0o - 0° as reference, one expects that the PA 
changes only marginally (less than 1°) within 10 yrs. A much 
more significant change of 10° (20°) would be expected after 
36 yrs (50 yrs), which would be detectable with future VLBI 
observations. 



limited (e.g. lCliffe et alJll996l) . we performed a new hydrody- 
namic simulation. Besides the general morphology, we aim for 
comparing the position of the newly discovered fainter stars 
with the column density variations caused by a precessing jet, 
which was beyond the scope of earlier studies. 

We use the version o f the ZEUS-3D code as modified 
by ISmith & RosenI ( l2003l) . which includes some molecular 
cooling and chemistry, as well as the ability to follow the 
molecular (H2) fraction. The large precession envisioned for 
the flow associated with NGC 7538 IRSl requires that the sim- 
ulation be performed on a very wide computational grid. Due to 
computational limits, we were restrained to use for this simula- 
tion a 3D Cartesian grid of 275 zones in each direction, where 
each zone spans 2 x 10'"* cm in each direction. This grid bal- 
ances the desire for some spatial resolution of the flow with 
the ability to simulate a sufficiently large part of the observed 
flow associated with NGC 7538 IRSl. Still, the total grid size 
(~ 0.018 pc) is smaller than the projected distance between IRS 
1 and A" -band feature C. 

Owing to the rather small physical size of the grid, we have 
chosen a nominal speed of 150 km s ', reduced from the in- 
ferred value of 250 km s ' for this source. The flow is pro- 
cessed with a nearly 30° precession angle, with the amplitude 
of the radial components of the velocity 0.55 of that of the axial 
component. The precession has a period of 120 years, which 
leads to 1.25 cycles during a grid crossing time. The flow is 
also pulsed, with a 30% amplitude and a 30 year period. This 
short period assists in the reproduction of the multiple knots 
of /T'-band emission near NGC 7538 IRSl. The jet flow also 
is sheared at the inlet, with the velocity at the jet radius 0.7 
that of the jet center. We have chosen a jet number density of 
10^ hydrogenic nuclei cm"^, while the ambient density is 10"^. 
The simulated jet radius is 4.0 x lO'^ cm (20 zones). Thus, the 
time-averaged mass flux is 2.6 x 10"^ Mq yr"', which is three 
orders of ma gnitude lower than the value as determined for the 
CO outflow JDavis et al.ll998l) . Similar calculations of the mo- 
mentum flux and kinetic energy flux, or mechanical luminosity, 
yield values of 3.8 x lO""* km s"' Mq yr"' and 4.7 Lq, respec- 
tively. 



4.4.4. Numerical molecular hydrodynamic 
simulations 

A large number of studies about the structure and evolu- 
tion of p recessing prqtostell a r jets can be found in litera- 
ture (e.g. iRaga & Birdll993t IVolker et al.lll999t iRaga et all 



w 

m 



l2004tlRosen & Sm ith-20041: ISmith & RoserJl2005al) . although 
most of these studies focus on jets from low-mass stars with 
rather narrow precession angles. As the number of simula- 
tions carried out for wide precession angles is much more 

* Positive values of 1// indicate an inclination out of the plane of the 
sky towards the observer. 

' For the sense of rotation, we follow the convention that counter- 
clockwise rotation (as measured from the source along the precession 
axis) corresponds to a positive sign of the phase (p. 



After the convolution with a PSF resembling the resolution 
obtained in real observations the H2 emission in our simula- 
tion shows a morphology which is similar to the one seen in 
the A" -band speckle image. In particular, the simulations might 
also explain features D and E as associated with the proposed 
precessing jet (compare Figure|8ji with Figure|3}. The simula- 
tions also show that the CO emission, which can be expected 
for a precessing jet at larger distances from the driving source, 
appears very smooth, which is also in accordance with the CO 
observations made for NGC 7538 IRSl. 



We note that the fainter stars e to n reported in Sec. 13.1.31 
are located in the region where the column density in our pre- 
cessing jet simulation appears particularly low (see Figure |8j^, 
left column), supporting the scenario proposed in Sec. 14.4.21 
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4.4.5. Possible precession meciianisms 

Several mechanisms have been proposed which can cause 
jet bending or jet precession, although most of them were 
established for low- and intermediate st ars and can cause 
prece s sion angles of only a few degrees ( iFendt & Z innecken 
Il998t lEisloffel & Mundtl Il997l) . For the cas e of h igh-mass 



stars a nd larger precession angles (0 ~ 40°), IShepherd et alJ 
( I2OO0I) summarized the three most promising concepts that 
could induce precession into circumstellar disks. We discuss 
how well these mechanisms can explain the observations of 
IRS 1 . For all cases, it is assumed that the outflow is launched 
close to the center of the disk and that a precession of the 
inner parts of t he disk will trans late into a precession of the 
collimated flow dBate et al.l200Cll) . 



1. Radiative-induced warping: lArmitage & Pring lel J 19971) 
suggested that geometrically thin, optically thick accretion 
disks can become unstable to warping if the incident radiation 
from the stellar source is strong enough. As this warping 
instability is expected to occur only at disk radii larger than a 
critical radius R^^i, we can estimate whether radiative-induced 
disk warping is expected at the inner part of the IRS l disk. 
Using a steflar mass of M ~ 30 Mq (Pestal ozzi et al.lE004ii 
a mass accretion rate of the order o f the mass outflo w rate 
Mace ~ Moutflow ~ 5.4 X 10"^^ MfTi yr-' ("Davis et al.'l998'), and a 
luminosity L j 9.6 x 10"^ Lcd ( Aka bane & Kuno 2005), we use 
equation 5 bvlArm itage & Pringld 1119971) and the assumptions 
listed in IShepherd et alJ 1I2OOOI) and obtain a critical radius 
^ciit ~ 200 pc. Since this is far beyond the inner edge of the 
disk where the jet collimation is expected to happen, it is 
very unlikely that the radiation emitted by the star or due to 
accretion causes any noticeable warping within the disk. 

2. Anisotropic accretion events: The impact/merging of 
(low mass) condensations can change the orientation of the 
disk angular momentum vector In such a dramatic event, 
angular momentum can be transferred from the impactor 
onto the accretion disk, potentially resulting in a net torque 
in the rotation of the disk. To estimate the precession angle, 
which could result from anisotropic accretion, very detailed 
assumptions about the disk, the impacting condensation, and 
their kinematics must be made. Since no data is available to 
est imate these quantities , we refer to the example computed 
bv lShepherd et al.l ( I2OO0I) and note that in extreme cases, such 
an accretion event could cause a sufficiently large precession 
angle in the case of NGC 7538 IRSl as well. However, in this 
scenario one would expect rather sudden changes in the jet 
direction rather than a smooth precession. 

3. Tidal interactions with a companion: Warping and pre- 
cession of the disk could be caused by tidal interactions with 
one or more companions on non-coplanar orbits. We assume 
the simplest case of a binary: with stellar masses Mp (primary) 
and Ms (secondary), an orbit with inclination / with respect to 
the disk plane, and a semimajor axis a. The mass ratio shall 
be denoted q = Mp/Mg and will be assumed as unity. Our 
observations place an upper limit on the separation of such a 
companion (see Figure|3i. Two cases can be considered: 



3a) circumprimary disk (a > ro): Because tidal torques 
woul d truncate the disk at ab out 0.3 times the binary separa- 
tion (iLubow & OgilvidbOOOl) . we obtain a lower limit for the 
binary separation (for a circular orbit), namely, a > 2 500 AU. 
However, a binary with such a large separation would be not 
suited to explain the observations since the orbital period 
would be > 2 X lO"* yrs (Mp + M , - 30 Mg^), imp lying a disk 
precession rate of > 4 x 10^ yrs (IB ate et al.ll200fll) . Assuming 
an extreme eccentricity might yield a short precession period 
of the order of 10^ yrs but implies strong, periodic interactions 
between the companion and the disk during each perihelion 
passage. As this would quickly distort and truncate the disk, we 
see this assumption contradicts the methanol maser structure, 
which suggests a smooth extension of the methanol layer from 
~ 290 AU to ~ 750 AU. 

3b) circumbinary disk (a < r,): Smoothed part ic le hy- 
drodynamic simulations by iLarwood & Papaloizoul (Il997h 
showed that a binary on a non-coplanar orbit with large 
inclination / could cause strong quasi-rigid body precession 
of the circumbinary disk (for ^ > 10) and strong warping, 
especially on the inner edge of the disk [q ~ 1). The same 
authors report that the disk precession frequency Wprec should 
be lower than the orbital frequency of the binary ^binary 
To make an order-of-magnitude estimation for the orbital 
period that would be expected for this hypothe tical IRSl 
binary system, one can assume Wpi-ec ~ Wbinai7/20 dBate et alJ 
I2OO0I) to obtain Pbinary ~ 14 yrs for the binary period, cor- 
responding to a separation of flbmaiy ~ 19 AU (~ 7 mas). 
This binary separation then puts a lower limit on the radius 
of the inner edge of the circumstellar disk. As this scenario 
can trigger the fast disk precession without truncating the 
extended disk structure traced by the methanol masers, we 
consider a circumbinary disk as the most plausible explanation. 



4.5. The IRS2 companion and flow interaction witii 
thelRS2UCH\\ region 

The spectral type o f IRS 2 was estimated to be 
04.5 dAkabane & Kunol l2005l) . corresponding to a lumi- 
nosity of ~ 6.4 X 10^ L0. Using the measured /iT'-band flux 
ratio, one can make rough estimates for the spectral type of the 
two components reported in Sec. 13.1.21 By assuming the total 
luminosity is attributed only to the two components, we obtain 
a spectral type of 05 for IRS 2a and 09 for IRS2fe (using the 
OB star luminosities from lVacca et alJl996l) . 

Within our images, the wide-opening angle outflow cone 
from IRSl seems to extend well out to IRS2. This offers an 
explanation for the s hock tracer line emit ting region that was 
imaged around IRS2 dBloomer et al.l 19981 see Figure|6l). The 
bowshock-like morphology of the [Fe II] and H2 emission sug- 
gests that the shock is excited from the south (which is roughly 
the direction towards IRSl). In the direction opposite IRSl, 
the [Fe II] and H2 emission even shows a cavity -like structure, 
which also appears in the 6 cm-radio continuum. lBloomer et alJ 
dl998l) suggested a stellar wind bowshock scenario, in which 
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IRS2 moves with a speed of ~10 km s ' towards the south- 
west through the ambient molecular cloud. We note that the 
morphology could also be explained by interaction between the 
IRSl outflow and IRS2 outflows. Based on its young age, 1RS2 
might also launch a powerful wind itself, causing the distinct 
shock zone which appears within the shock tracer emission (see 
Figure |6j) and which is also detectable in our /T'-band image 
(arc-like morphology between features G and H). 

4.6. Outflow Structures from IRS1 at larger spatial 
scales 

Figure |6j; to c/ shows mosaics of the vicinity of IRS 1 in the 
four IRAC bands. Although IRS 1 and 1RS2 appear saturated in 
these images (shown in logarithmic scaling) and banding (verti- 
cal and horizontal stripes produced by IRS 1 and IRS9) appears 
especially in the 5.8 and 8.0 yum images, structures potentially 
related to IRS 1 can be observed. About 40" towards the south- 
east of IRSl, a bowshock structure can be se en, which is also 
present in the H2 image bv lPavis et al.l (119981) . This bowshock 
points in the same direction as the redshifted lobe of the CO 
outflow (see Figure|6f) and just opposite to the outflow direc- 
tion identified in our speckle image at small scales. Thus, it is 
possible that this bow traces the southeastern part of the IRS 1 
outflow. In our speckle image, the inner part of this southeast- 
ern outflow is not visible; likely a result of strong intervening 
extinction. 

Furthermore, it is interesting to note that the "ridge" con- 
necting the IRS 1-3 cluster with IRS4 and IRS5 follows the 
western wall of the outflow direction identified in our speckle 
image (see top of Figure |6}. It is possible that the total extent 
of the IRSl outflow also reaches much further northwest than 
the structure seen in the speckle image, contributing to the ex- 
citation of the western part of the bubble seen in the IRAC 
bands and the shocks in the H2 image by IPavis et al.l lll998i 
see Figure|6l2). 

5. Evidence for triggered star formation in thie 
NGC 7538 star forming region 

It has been proposed by many authors that star formation 
seems to p ropagate southeastwa rd s thr oughou t the NGC 7538 
complex JWerner et all Il979t JMcCaughre an et al.l Il99lt 
lOjha et al.ll2004l) . This is indicated by the spatial arrangement 
of the members of this star forming region, which also seems 
to agree with the expected evolutionary sequence: starting 
about 3' northwest of IRSl, O stars located in the H II region 
represent the most developed evolutionary state, followed by 
the IRS 1-3 cluster and their associated UCHII regions, with 
the compact reflection nebula around IRS9 representing the 
youngest membe r of this star form ation site. In agreement 
with this picture, Balog et al.l (l2004l) measured the reddening 
of stars throughout NGC 7538 and found a gradient in red- 
dening with the most heavily reddened sources in the southeast. 

The presented IRAC images can also be interpreted in 
support of this scenario since the "ridge"-like feature connect- 
ing IRS 1-3, IRS4, IRS5 seems to trace the interface between 



the northeastern bubble (visible at NIR/MIR wavelengths) 
and the submillimeter bubble, which appear s in the 450 //m 
and 850 //m-maps by iReid & Wilsonl ( l2005l see Figure |2}. 
This suggests that in NGC 7538, star formation was triggered 
by the compression of gas just at the interface layer of these 
expanding bubbles, sequentially initiating the formation of the 
observed chain of infrared sources. 



lOiha et al.l ll2004l) suggested that IRS6, the most luminous 
source in the NGC 7538 region, might be the main excit- 
ing source responsible for the optical H II region. Inspecting 
the IRAC color composites, this scenario is supported by the 
morphology of the bright, curved structure west of IRS6. In 
the IRAC 8 //m band (red in Figure I2J7), this structure ap- 
pears particularly prominent. As it is known that emission 
in this IRAC band is often associated with PAHs, this sug- 
gests that this region is illuminated by strong UV radia- 
tion from IRS6. Other features, such as the conical structure 
around 2MASS 23 1 35808H-6 1 30484 and the structure northeast 
of IRS7, also show a symmetry towards IRS6. 

6. Summary and Conclusions 

Bispectrum speckle interferometry and archival Spitzer/IRAC 
imaging of the massive protostars NGC 7538 IRS 1/2 and their 
vicinity are presented. We summarize our results as follows: 

1. The clumpy, fan-shaped structure seen in our speckle im- 
ages most likely traces recent outflow activity from IRSl, 
consistent with the direction of the blue-shifted lobe of the 
known CO outflow. A bowshock structure noticeable in 
the IRAC images ~ 40" southeast of IRS 1 suggests that 
the total extent of the outflow might be several parsecs. 
The outflow might have also contributed to shaping and 
exciting the bubble-like structure, which is prominent in 
all four IRAC bands (although contributions from several 
other sources, especially IRS6, are also evident). 

2. A companion around the high-mass star NGC 7538 IRS2 
was discovered. Furthermore, we see indications for inter- 
actions between the IRS 1 flow and outflows or stellar winds 
from IRS2 (nebulosity surrounding IRS2). 

3. A jet precession model seems suitable to describe the fea- 
tures observed within our NIR images, simultaneously ex- 
plaining the misalignment between the putative methanol 
maser disk, the UCH II region, and the outflow tracers de- 
tected at larger scales (CO). A simple analytic precession 
model was used to extract order-of-magnitude estimates for 
the precession parameters. Using these we estimate tidal in- 
teraction of a close binary system with a circumbinary disk 
as the most plausible gyroscopic mechanism, which is trig- 
gering the precession. 

4. The presented molecular hydrodynamic simulations can re- 
produce some of the fine-structure observed in our NIR 
images and indicate that the arrangement of the detected 
fainter stars might be explained as a column-density effect, 
caused by the proposed precessing jet. 

5. The prominent sites of ongoing high-mass star formation 
in NGC 7538 seem to be located just at the interface be- 
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tween two bubble-like structures — one is visible in the 
presented IRAC images, the other traced by submillimeter 
observations. The gas compression caused by the expan- 
sion of these bubbles might have triggered star formation 
in this region. 

While it is well estabUshed that the outflows of HMPOs 
generally appear less collimated than those of their low-mass 
counterparts, the recent discovery of evidence for outflow 
precession for an increasing number of massive YSOs might 
indicate a common launching mechanism for all outflow 
driving sources of all stellar masses. The observed w idening in 
HMP O outflows might be due to selection effects dShepherdl 
I2OO5I) and/or precession of a collimated jet. T he large preces- 
sion angles reported for IR AS 20 1 26-h4 1 04 dShepherd et"ai] 
l2000l) . S1 40 IRSl llWeigelt et al.1 l2002h . IRAS 
23151H-5912 JWeigelt et al.ll2005h . and now NGC 7538 IRSl 
(this paper) might point towards a rather dramatic precession 
mechanism, maybe the presence of very close, high-mass 
companions on non-coplanar orbits. 

We strongly encourage further observations of IRSl, es- 
pecially to detect potential companions either by near-infrared 
long-baseline interferometry or radial velocity measurements. 
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Fig. 2. Figure a) shows the Spitzer/IRAC 4.5 /im image with the position of the infrared sources (IRSl to 11 and 
2MASS sources) a nd H7O masers marked. Furthermore, the position of the sub-millimeter (450 yum and 850 /im) 
clumps reported by iReid & Wilsonl (l2005l) are shown. The position of the 2MASS sources 2MASS 23 135808H-6 130484, 
2MASS23131660H-6128017, 2MASS23134351H-6129372, 2MASS23144651H-6129397, 2MASS23131691H-6129076, 
2MASS23130929H-6128184, and 2MASS23133184H-6125161 are labeled explicitly. Figures h) and c) show color-composites 
produced with two triplets of the four IRAC bands. The intensity of each image was scaled logarithmically. North is up and east 
to the left. 
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Fig. 3. Bispectrum speckle image with identified point sources (triangles) atop marked. The astrometry for the point-sources was 
performed using the high-resolution BTA image, whereas the image shown was reconstructed from MMT data. The contours 
trace 0.25%, 0.5%, 0.75%, 1.0%, 1.25%, and 1.5% of the peak intensity. The inset on the upper left shows a reconstruction of 
the vicinity of IRS2 using a resolution of 80 mas (BTA data). In the lower left, IRSl is shown using a different color table, 
emphasizing the elongation of the IRSl Airy disk (MMT data) overplotted with the 15 GHz radio continuum (the contours 
sho w -1, 1, 2.5, 5, 10, 20, ...,90% o f the peak flux) and the position of the OH (circles) and methanol (crosses) masers (image 
from lHutawarakorn & CoherJ2003l u sing data from | naume et alJll995l) . In the lower right we show the integrated brightness of 
the methanol masers as presented bv lPestalozzi et alJ (120041 



contour levels of 1, 3, 5, 10, 30, 50, 70, and 90% of the peak flux 



density are shown). 
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Fig. 5. Illustration showing the 
outflow directions in the vari- 
ou s tracers. Th e CO contours 
by iKameva et alJ lll989t red and 
blue) are overl aid on the H? map 
(greyscale) by IPavis et alJ ( Il998l) . 
The orientation of the conjectural 
methanol maser disk (green), the 
fan-shaped structure detected in 
our A" -band image (orange), and 
the averaged direction of H2 (red 
arcs) are shown schematically. 
The arrows indicate the direction 
prependicular to the alignment of 
the methanol masers (green), the 
orientation of the inner (< 0.5") 
and outer (> 0.5") core detected 
in the 1.0, 1.3, and 6.0 cm radio 
continuum (white), and the direc- 
tion along which the IRSl Airy 
disk w as found to be elongated 
(MIR: be Buizer & Minied l2005t 
NIR: this paper). 
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Fig. 6. Mosaic showing the IRS 1-3 cluster at various w avelengths. Bes ide the speckle /f'-band im age (al so marked as red box) 

^li COl E^evaetalJliosi COVlRloomer etTilflQoi J, H, K, H2, [Fell]), 

lDavisetal.l(ll998l Uo). and lPe Buizer & Minieil (l2Q05[ 1 1.7 //m, 18.3 y^m) was incorporated. 



and IRAC irnages, data fron ilScoville et alJ 1119861 COuKameva et i 



Kraus et al.: Outflows from the high-mass protostars NGC 7538 IRSl/2 



19 










ff^^ 



►•West 



2 - 



1 - 



1 


1 1 1 1 1 1 

_ Precession Axis 

" ^ Radio ContrnuLrm 


1 






— 


/ Jet Trajeetofv 




1 


^^\^~^ 


• 


•t 




c 


■ 


^k^ 


' 


a 


i' 




^ 


w 

1 1 1 1 1 t 


1 



■2 -3 

RA offset ("] 



-5 



Fig. 7. Left: Illustration of the analytic precession model presented in Sec. 14.4.31 Right: MMT speckle image overplotted with 
the trajectory of ejecta from a precessing outflow projected onto the plane of the sky (thick blue line) as described by the analytic 
precession model. For the counter-clockwise precession, the parameters f prec = 280 yrs, = 40°, a - -30°, i/' = 0° (precession 
axis within the plane of the sky), and 0o - 0° were used. In order to simulate a finite collimation of the flow, we varied both 
a and i/r by +5°, yielding the traj ectory given by the thin blue line s. The red contours show the 15 GHz radio continuum map 
bv lHutawarakorn & Cohenl2003l ("using data bv lOaume et alJl995l) . 
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Fig. 8. flj, fej: Synthetic images (H2 1 ^ 0, 2 -^ 1; CO R(l), R(5)) from our ZEUS-3D molecular hydrodynamic simulation, 
shown for two different projections. In Figure a) the images are also shown after a convolution with a Gaussian which roughly 
resembles the resolution obtained in our speckle observation (Figure|3}- Figure c) shows the total gas column density and the H2 
column density for the same projections. Finally, in d) we show channel maps of the CO outflow for four velocity bins. Each 
velocity bin has a width of 5 km s"' and the central velocity is given by the number on the left of each image (in km s"'). The 
two numbers on the right of each image indicate the log of the maximum integrated luminosity in any single element in the image 
{Top, in erg s ') and the log of the total integrated luminosity in the entire velocity bin {Bottom). 

The angle (f> gives the angle between the z axis and the LOS, corresponding to a rotation around the precession axis {x axis). 
Please see the text (Sec. 14.4.41 for a description of the complete model parameters. 



